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TECHNICAL NOTE 2408 

APPROXIMATE DESIGN METHOD FOR HIGH- SOLIDITY BLADE ELEMENTS IN 

COMPRESSORS AND TURBINES 


By John D. Stanitz 


SUMMARY 

An approximate blade-element design method is developed for com- 
pressible or incompressible nonviscous flow in high-solidity stators or 
rotors of axial-, radial-, or mixed-flow compressors, turb ine s, or two- 
dimensional cascades. The method is based upon channel- type flow 
between blade elements on a specified surface of revolution that lies 
between the hub and shroud (casing) and is concentric with the axis of 
the compressor or turbine. The blade element is designed for prescribed 
velocities along the blade-element profile as a function of distance 
along meridional lines on the surface of revolution. The method is 
limited, because of assumptions, to prescribed velocities that result in 
blade-element profiles with gradual variations in thickness and with 
mi . n1 . Tnu m radii of curvature at least approximately equal to the channel 
width between profiles. 

Two numerical examples are presented: The first example is the 

design of a blade-element profile for a plane two-dimensional cascade in 
compressible flow with prescribed velocities along the profile; the 
second example is the design of a blade element for the impeller of a 
mixed-flow centrifugal compressor. In both examples the design method 
has been checked by comparing the prescribed velocity distribution with 
the velocity distribution obtained by stream?- filament methods from the 
resulting blade-element design. 


INTRODUCTION 

Most aerodynamic losses in compressors and turbines result from 
the viscosity and compressibility of the fluid. (Other losses, for 
example, are associated with the trailing vortex sheet that sheds from 
the trailing edge of the blade if the circulation around the blade 
elements varies along the span of the blade.) Because the fluid is 
viscous, losses result from boundary- layer friction and separation and 
from secondary flows associated with the boundary layer. The magnitude 
of these losses depends on the boundary- layer thic kne ss, the size of 
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which, is controlled by the velocity distribution just outside the 
boundary layer. In particular, if the velocity decelerates too rapidly,, 
the boundary layer thickens and separates, causing large mi xl ng losses. 
Also, because the fluid is compressible, significant shock losses result 
if the maximum relative velocity on the blade surfaces appreciably 
exceeds the local speed of sound. In addition, if the relative velocity 
exceeds the local speed of sound for too great a distance across the 
channel between blades, "choke-flow" conditions develop and the blade 
row cannot pass the desired quantity of fluid. By proper selection of 
the velocity distribution relative to the blade surfaces, secondary- flow 
losses and boundary- layer friction and separation losses can be reduced, 
and shock losses and choke-flow conditions can be e liminat ed. With the 
assumption that desirable velocity distributions can be determined from 
boundary- layer theory and so forth, the problem of efficient compressor 
and turbine design is resolved into the development of design methods 
for blades that have prescribed velocity distributions. 

In the past, a number of design methods have been developed for 
plane, two-dimensional cascades with prescribed velocities along the 
blade surfaces for incompressible flow (references 1 to 3, for example) 
and for compressible flow (references 4 and 5) . These cascades are 
limited in application to axial- flow compressors and turbines. A method 
(reference 6) has also been developed for the design of radial-flow 
centrifugal impellers for incompressible flow and for a prescribed 
variation in work input with radius . Recently, work has been done 
(references 7 and 8, for example) on the general problem of blade 
design for mixed- flow compressors or turbines in which various inde- 
pendent variables are prescribed, but usually the relative velocity along 
the blade surfaces is not. 

The work presented herein is concerned with an approximate design 
method developed at the HACA Lewis laboratory for blade elements with 
prescribed velocities along the profile . The method is developed for 
elements of high-solidity blades in compressors and turbines or for 
plane, two-dimensional cascades. The blade-element profile lies on a 
specified surface of revolution that is concentric with the axis of 
the compressor or turbine and lies on or between the hub and shroud 
(casing) . The blade-element-design method is for compressible, non- 
viscous fluids and, because of assumptions required by the method, is 
limited to prescribed velocities that result in blade- element profiles 
with gradual variations in thickness and with miinimum radii of curvature 
at least approximately equal to the cha nn el width between profiles. 

(The nose and tail regions of the blade-element profile are excluded 
from the approximate blade- element- design method of this report.) 
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THEORT OF DESIGN ME2EHQD 

She approximate blade- element- design method is developed in two 
parts. The first part, called the first approximation, is simple and 
rapid and appears to he satisfactory for very high-solidity blades, 
such as exist in radial- and mixed- flow centrifugal compressors. The 
second part, called the second approximation, improves on the assump- 
tions of the first part and appears to be satisfactory for blade 
solidities such as exist in present-day turbine nozzles. 


Pre li minary Considerations 

Analysis problem. - In the aerodynamic theory of compressors and 
turbines there are two types of problem: first, the direct or 

analysis problem in which the velocities are deter mine d for a given 
blade shape; and second, the inverse or design problem in which the 
blade shape is determined for prescribed velocities along the blade 
surfaces. It is the inverse or design problem with which this report 
is concerned, but by way of introduction to the design method the 
direct or analysis problem will be discussed first. 

One possible means of analyzing the flow through a blade row 
(stator or rotor) with specified geometry is as follows: Consider the 

flow of an ideal, compressible fluid through an arbitrary channel 
between blades, such as shown in figure 1. The fluid is free to follow 
whatever path the pressure and inertia forces require of it. If 7 how- 
ever, the number of blades in the blade row approaches infinity, the 
space between blades approaches zero and the path of the fluid is 
restricted to the curved, mean surface of the blade. (The blades 
become very thin so that the two surfaces of each blade approach a 
mean surface.) Under this assumption of axial symmetry (references 7 
and 9, for example) the fluid flow is reduced from three-dimensional 
motion in the passage between blades with finite spacing to two- 
dimensional motion on the curved, mean blade surface. The streamlines 
of this two-dimensional motion can be projected on the meridional plane 
(axial-radial plane) as shown in figure 2. Ruden (reference 10) has 
shown that, provided the blades are not too widely spaced, axial- 
symmetry solutions give a good picture of the mean flow between blades. 

For finite blade spacing, flow conditions vary from blade to blade 
( circumferentially about the axis of the compressor or turbine) as well 
as from hub to shroud. In order to investigate the variation from 
blade to blade, it is assumed that the motion of any fluid particle 
bounded by adjacent streamlines in the meridional plane (fig. 2) is 
restricted to the annulus generated by rotating these adjacent 
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streamlines about the axis of the compressor or turhiue. If the adjacent 
streamlines are sufficiently close together, flow co ndi tions in the 
annulu s can he considered uniform normal to the mean surface of revo- 
lution -wi thin the annulus. Thus the fluid motion is reduced to two- 
dimensional flow on the mean surface of revolution (fig. 3) generated 
by rotating the center line between the adjacent streamlines in the 
meridional plane (fig. 2) about the axis of the compressor or turbine. 
Blade- to-blade solutions of this type are given in references 11 and 12. 

Blade-to-blade solutions can be obtained for every mean surface of 
revolution generated by the center lines between adjacent streamlines in 
the meridional plane (fig. 2) . Therefore, flow conditions can be deter- 
mined approximately throughout the passage. The resulting quasi three- 
dimensional solution of the direct or analysis problem was obtained by 
the combination of two types of two-dimensional solution; axial- symmetry 
solutions in the meridional plane and blade-to-blade solutions on sur- 
faces of revolution. Such a combination of solutions prohibits the pos- 
sibility of a corkscrew path, which the fluid might follow In an exact 
three-dimensional solution, but it can be expected to give a better 
picture of the flow than does any two-dimensional solution alone. 

Design problem. - It is proposed to solve the inverse or design 
problem by an approach similar, but reversed, to that outlined for the 
direct or analysis problem in the previous section. First, the profile 
of the blade element on one surface of revolution (fig. 3) will be 
determined for prescribed velocities along the profile. In order to 
accomplish this first phase of the blade-design method, the Bhape of 
the surface of revolution (fig. 3) is specified together with the var- 
iation in height of the fluid particles (or blade element) along the 
surface. The specified surface of revolution Is generated by a 
prescribed curve that is the center line between two adjacent stream- 
lines in the meridional plane. (The curve should be smooth with large 
radii of curvature, but Is otherwise assumed arbitrary in shape.) The 
specified variation in height of the fluid particles (or blade element) 
along the surface of revolution is given by the prescribed spacing of 
the adjacent meridional streamlines that generate the annulus around 
the surface of revolution. (The variation in meridional streamline 
spacing along the center line should be smooth and slowly varying but 
is otherwise assumed arbitrary.) 

After the profile of the blade element on the surface of revo- 
lution has been determined by design methods that will be developed 
in this report, the shape of the blade surfaces that cut through 
the profile of the blade element on the surface of revolution must 
be obtained and the hub. and shroud contours of the blade in the 
meridional plane must be determined. A suggested method to accom- 
plish this second phase of the design method (not the subject of this 
report) is first to prescribe the blade surfaces. (These surfaces 
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pass through the blade- element profile on the specified surface of revo- 
lution and extend to the hub and shroud casings, "which are as yet 
■undetermined. The prescribed shape of the blade surfaces determines 
the blade taper and the thickness distribution everywhere except on 
the specified surface of revolution.) The hub and shroud profiles are 
then determined by a solution to the inverse or design problem in the 
meridional (axial- radial) plane starting from the specified shape of the 
center line between adjacent streamlines in the meridional plane and 
the specified spacing between the adjacent streamlines and solving for 
the variation in meridional stream function outward and inward along 
normals to the meridional streamlines until the meridional strea mli nes 
along the root and tip (hub and shroud) of the blade are determined. 

The specified spacing of the adjacent stream l i n es together with 
the known blade-element-profile thickness distribution determines the 
variation in meridional velocity along the specified center line between 
the streamlines. This arbitrary specification of the magnitude of the 
velocity along a specified stream path (the center line) overdetermines 
the problem and leads to singularities in the flow field. If, however, 
the prescribed velocity is reasonably consistent with the prescribed 
shape of the center line and if the velocity is prescribed along the 
center line at a finite number of points only, the hub and shroud 
(casing) can be determined (second phase of blade-design method) by 
finite-difference methods that in general are not troubled by the 
singularities. The resulting hub and shroud contours can be checked by 
a direct solution for the shape of the streamline corresponding to the 
prescribed center line and for the spacing of the adjacent streamlines. 
If the agreement with the prescribed center line and spacing of the 
adjacent streamlines is not good, either a new center line or a new 
meridional velocity distribution (spacing of the adjacent streamlines) 
or both are prescribed and the entire blade-design method is repeated 
until satisfactory agreement is achieved. 

Thus the entire blade geometry can be determined for specified 
velocities along the profile of a blade element on one surface of 
revolution, which surface is specified together with the arbitrary (but 
smooth) variation in height of the fluid particles (or blade element) 
along the surface. Only the first phase of the proposed blade-design 
method will be considered in this report; that is, only the profile 
of a blade element on one surface of revolution will be determined for 
prescribed velocities along the blade profile. If the hub-shroud 
(root- tip) ratio of the blade is sufficiently large and if the speci- 
fied shape of the center line between meridional streamlines has suffi- 
ciently small curvature, the second phase of the design method can be 
neglected, as is customarily done in the application of plane, two- 
dimensional cascade data to the design of axial-flow compressors and 
turbines, for example. 
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Coordinates . - The cylindrical coord in ates R, 0, and Z are shown 
in figure 3. (All symbols are defined in appendix A.) These coordinates 
are dimensionless, the linear coordinates R and Z having "been divided 
hy a characteristic length equal to the chord c of the "blade element 
(the exact definition of c will he given later) . The coordinate system 
is oriented with its Z-axI s along the axis of -the compressor or turbine. 
The coordinates are fixed relative to the blades so that the coordinate 
system is absolute if the blades are stationary (stators) and relative 
if the blades are rotating (rotors) . The angular velocity (a for 
rotors is always positive (right-hand rule) about the Z-axis , as shown 
in figure 3. 

An elemental distance dS in the direction of flow (that is, 
coinciding with the velocity vector) has components dR, Rd0, and dZ 
(fig. 3). The projection of dS on the meridional plane is given by 
dM in figure 3. The elemental distances dS and dM help to define 
two angles a and p, where, from figure 3, 


dR = dM sin a, 

(la) 

dZ = dM cos a 

(lb) 

and 


dM = dS cos P 

(2a) 

Rd0 = dS sin p 

(2b) 

Because dS and dM are always positive and finite, 
requires that 

equation (2a) 

2 <P< 2 


From equation (2b) and figure 3, p>0 for positive values of Rd0, and 
P< 0 for negative values of Rd0. Because both dS and dM lie on 
the surface of revolution, the angle p is the flow direction measured 
on this surface and is positive to the right when the exterior of the 
surface is viewed in the direction of dM (fig. 3) . 

In general, dZ can have both positive and negative values. In any 
blade row of most compressors and turbine?, however, dZ will always 
have the same sign, which in this report shall be considered positive 
(or zero) , so that from equation (lb) , 

-it- % 
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For inflow machines (fig. 4(a)), dR is negative in the direction of 
<3M and. for outflow machines (fig. 4(h)), dR is positive, so that 
from equation (la). 


— < a < 0 (inflow machine) 

0 < a < ^ (outflow ma chin e) 

For axial -flow machines, dR approaches 0 so that, from equation (la), 
a is approximately 0. 

Velocity components . - The velocity Q relative to the coordinate 
system has components Qp, %, and Q 2 in the R-, 0-, and Z- 

directions, respectively (fig. 3) . These velocities are dimensionless, 
having heen divided hy a characteristic velocity equal to the absolute 
stagnation speed of sound a 0 upstream of the blade row, where 


a Q 2 = rgRT 0 (3) 

in which R is the gas constant, r is the ratio of specifip heats, 
g is the acceleration due to gravity, T is the Btatic (stream) tem- 
perature, and the subscript o indicates stagnation conditions upstream 
of the blade row. For rotors, the tip speed of a blade element % is 
likewise dimensionless and is defined by 


0)RnC 

°T - -57- W 

Thus, if 

E * ■ 5 < 5 > 

the tangential velocity of the blade at any radius R is equal to R*tJrji 

and the absolute tangential velocity of the fluid is equal to 

(R*UT + Qq) . For stator blades, od is zero and therefore Uij is zero. 
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In addition to the three components of the relative velocity Or, 
Qq, and Q z , it is convenient to define a fourth velocity component, 
lying in the meridional plane and defined as the meridional velo- 
city Qy, "where from figure 3, 


Also from figure 3, 


and 


% = Q cos p 

(6a) 

% = % sin a 

(6b) 

Q# = Q sin 0 

(6c) 

Qg = 0^ cos a 

(6d) 


Fluid particle . - As shown in figure 5, a fluid particle on a 
surface of revolution has the dimensions RdS, dM, and AH. The dif- 
ferential angle d0 is measured in a plane normal to the Z-axis 
(fig. 5) and the differential distance dM is related to the differ- 
ential radius dR by 


dM = 


dR 

sin a 


(la) 


The incremental height AH (figs. 2 and 5) is measured normall. to the 
surface of revolution and is determined by the specified spacing of the 
adjacent streamlines in the meridional plane (fig. 2) . The incremental 
height AH is dimensionless, having been divided by the characteristic 
length c. The height ratio H* of the fluid particle (fig. 5) is 
defined as the ratio of the incremental height AH of the blade element 
at radius R to the incremental height (AH)rj of the blade element 
at Rrp. 


Thermodynamic relations . - From the general energy equation, the 
static (stream) temperature can be shown to be related to the relative 
velocity Q by (reference ll) 


^ = i + ^ [(E%) 2 - Q 2 “ 2% Nj] 


(7) 
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■where the subscript U refers to conditions upstream of the "blade row 
and 'where X is the whirl ratio (absolute moment of m omen tum divided 
by c%a Q ) given by 


X = R^lfUrji + Q 0 ) (8) 

The pressure ratio P and the density ratio p/p Q are likewise 
related to the relative velocity Q by 


0 ) 


( 10 ) 


First Approximation 

The first approximation gives a rapid means of computing the blade- 
element profile on any specified, surface of revolution (with specified 
variation in H*) for prescribed velocities along the profile. The 
method is based upon one-dimensional, channel -type flow between blade 
elements. 

Assumptions . - For the first approximation, it is assumed that on 
the surface of revolution in the direction of 6 (that is, for a given 
value of M) flow conditions may be considered uniform and equal to 
mass-weighted average values. This assumption is usual for one- 
dimensional, channel -type analyses, except that conditions have been 
considered uniform in the direction of 9 rather than along a line 
normal to the average direction of flow. In particular, it is assumed 
that at a given value of M: (l) The mass -weighted average flow direc- 
tion is equal to the direction of the blade-element camber line, and (2) 
the mass-weighted average velocity is equal to one-half the sum of the 
prescribed velocities at the two blade-element surfaces. Assumption (l) 
is approximately true for high- solidity blade elements with gradual 
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variations in profile thickness and is exactly true for infinite soli- 
dities. Ass ump tion (2) is approximately true if the velocity distri- 
bution across the channel in the 0-directian is linear, a condition 
that is approached for high-solidity blade -element profiles with mini- 
mum radii of curvature at least approximately equal to the channel 
width between profiles (see fig. 15 of reference 11 and figs. 7 and 16 
of reference 13) . 

Flow field. - A view on the SM-plane of the passage between two 
blade-element profiles is shown in figure 6. The flow field considered 
by the blade-element-design method lies between the profiles of the 
blade elements and between two parallel lines normal to the direction 
of Qm and located close to the leading and trailing edges of the 
blades, as shown in figure 6. The immediate nose and tail regions have 
been excluded from the flow field considered by the blade-element— design 
method, because the small radii of curvature along the blade-element 
profile in these regions result in distributions of velocity and flow 
direction circumferentially across the channel between blade elements 
tha t do not satisfy the assumptions of the design method in these 
regions. In the final design, the nose and tail contours are rounded 
off in a manner guided by experience. 

By definition, the characteristic length c (blade chord) is equal 
to the distance (not dimensionless) between the n- and t-boundaries 
(fig. 6) along a meridional line. Thus, the meridional distance M 
(dimensionless) when measured from the n-boundary is equal to 1.0 at 
the t -boundary. 

The blade solidity a is defined by 

= RpCAS) (11) 

which, from figure 6, is the ratio of the blade chord c to the blade 
spacing at the tip of the blade element cErj (A0 ) where (A0) is the 
angular spacing of the blades. 

The ratio of channel width H(0i-0o) to blade spacing E(A0) at 
a given value of M is defined by 


g i- e o 

A0 


( 12 ) 
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where the subscripts 0 and 1 refer to the blade surfaces at the left 
and right of the channel walls, respectively, when viewed in the direc- 
tion of Thus, from equations (5), (ll), and (l2), the blade 

spacing (fig. 6) at a given value of M is equal to R *fa and the 
cha nn el width is equal to 6R */cr. 

Prescribed velocities along blade-element contour . - In the pro- 
posed design method, the velocities Qq and along the blade- 

element profile on the surface of revolution are prescribed as func- 
tions of M. 

Outline of method . - The object of the design method is to deter- 
mine the blade-element profile (on a specified surface of revolution 
with specified variations in, H*) for prescribed velocities along the 
blade-element contour and for prescribed flow conditions (whirl ratio 
and flow rate) upstream and downstream of the blade-element row. An 
inspection of figure 6 shows that the blade-element profile is com- 
pletely described by the solidity a and by 


6 = 5(M) 

(13) 

P = P(M) 

(14) 


These quantities a, 5, and p can be determined (approximately) 
from the prescribed velocities Qq and Qq and from the prescribed 
upstream and downstream conditions in a manner outlined as follows: 

(a) The Bolidity cr is determined by the blade-element spacing, which, 
for a given flow rate, depends on the specified change in whirl ratio 
upstream and downstream of the blade-element row and on the total blade- 
element torque; (b) the flow direction p at each value of M is 
obtained by equation (6c) from the average values of Q and Qg. The 
average value of Q is by assumption equal to one-half the sum of Qo 
and Qq. The average value of Qg at M is determined by the blade- 
element torque between M = 0 and M; and (c) the ratio 6 is deter- 
mined from continuity considerations. 

Solidity a . - The blade-element solidity is determined frcm the 
condition that the rate of change of moment of momentum of the fluid 
flow between two blade elements must equal the torque of the blade ele- 
ments on the fluid. For a fluid strip of width cdM between blades 
(fig. 7) the torque of the blade elements is given by 
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■where 


- |p 0 (AP)j [e(AH) T E*| 




AP = P-l - 


P 


0 


(15) 


(16) 


In ■which P is determ in ed from equation (9) and from the specified 
velocities and Qq. 


The rate of change of moment of momentum of the fluid flow across 
the strip in figure 7 is given by 


•*¥>0 f| ^ar « ' (17) 

■where h is the number of "blade elements (or passages "between "blade 

elements) , ^ is the incremental flow rate across the fluid strip 

(fig. 7) so that Aw is the total incremental flow rate along the sur- 
face of revolution, and the subscript av refers to the mass -weighted 
average value in the circumferential direction across the channel 
"between "blade elements at a given value of M. 

Equating the quantities given "by (l5) and (17) and introducing the 
perfect gas law (p = pBT) give 


d *av 

dM 

where the flow coefficient 


= — H*R* AP 

rep 

is defined "by 


(18) 


9 


1 Aw 

p a (Af ) m 
K o o v T 


(19) 


where (Af)rp, the total ("blocked plus unblocked) incremental annulus 
area normal to the direction of ^ at the blade-element tip, is given 
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(Af) T = [b(A0]] [cr£] [c(AH)^] (20) 

Tims Aw/(Af)rji in equation (19) is the flow rate per unit total area 
at the "blade -element tip. 

The whirl ratio X aT varies from ^ to ^ as M varies from 
0 to 1.0 so that equation (l8) integrates to give 


a = 


- r qffip ~ X u) 

B*R* (AP) dM 


( 21 ) 


Jft * 

where H , R , and AP are known functions of M so that cr is 
directly determined hy equation (21). 

Average relative tangential velocity ( Qg) av « - The mass -weighted 
aver ge relative tangential velocity at any specified value of M is 
c , ' 4 -sined hy integrating equation (l8) from M = 0 to M so that 


“ \j yep 


f\M 


H*R* (AP) dM 


°0 

lich, from equation (8), becomes 


(Qe) 


_ x u ^ a 

av r* " R TJ T - yqe* 


nM 


E*R* (AP) dM 


UO 


( 22 ) 


where n is known from equation (21) so that (Qo) a v i® directly 
determined hy equation (22) . 

Flow, direction ft . - From equation (6c) and from the assumption 
that Q and p may he considered constant for any given value of M 
in the direction of 0 between blade elements on the surface of revolu- 
tion. 
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P = sin"l 


(Qf^av 

%.v 


(23) 


■where, as all ready stated -under the assumptions, the average value 
of Q is 


90 * <*L 
2 

Equations (22), (23), and (24) determine the flow direction P as a 
function of M. 

Channel-width ratio 8 . - The channel -width ratio is determined by 
cont inui ty considerations. The incremental flow rate Aw/h across the 
fluid strip in figure 7 is given hy 



T = [ oB ^i- e o] [ e ^) T ^ [ p ° t“] M oos “ 

from which, together with equations (l2), (l9), and (20), 

8 - j- 2 (25) 

bV 22 cos 8 

The ratio p av /p 0 is given hy equation (10) with Q equal to Q^. 

Equation (25) dete rmin es 6 as a function of M. Upstream and down- 
stream of the blade-element row, 8 equals 1 so that from equation (25), 


cp = (r*H* ■£- Q cos pL = (r*H* f - Q cos pV 


( 26 ) 


which determines 

H*R* Q cos P 
^o 


a relation that must he satisfied between 
upstream and downstream of the blade-element row. 
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Blade -element -profile coordinates . - The value of a has "been 
determined, and 5 and p are known functions of M so that the 
design of the "blade-element profile on the surface of revolution is 
complete (see equations (l3) and (l4)). The coor dina tes of the "blade- 
element profile on the surface of revolution are obt ain ed "by (l) com- 
puting the shape of the profile camber line (fig. 6) on the surface of 
revolution, and- (2) distributing the blade thickness equally on either 
side of the camber line. 

The shape of the profile camber line is given by R and 0 as 
functions of M. Because R (or R*) is a known function of M, it 
only remains to determine R0 as a function of M in order to obtain 
the profile camber line. Along the camber line. 


d(R9) _ d0 dR 

e m. 

which becomes, when combined with equations (la), (2a), and (2b), 


= tan P + 0 sin a (27) 

The change in R0 from 0 at M = 0 to R0 at M is obtained by 
integrating equation (27) to yield 


rU 

R0 = (tan p + 9 sin a) <3M (28) 

Jo 

Because R (or R*) is a known (specified) function of M, equa- 
tion (28) determines the blade-element camber line on the surface of 
revolution. 

The blade thickness in the 0-iLrection is obtained from the 
channel -width ratio 8 and is given by 



(29) 


R(A0) (1-8) 
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[phis expression, for the "blade thickness is dimensionless, and is 
expressed, as usual, in units of characteristic length c. The blade 
thickness is equally distributed cn either side of the camber line. 

The hlade-element spacing on the surface of revolution is given 
directly by the solidity cr. 


Second Approximation 

The second approximation modifies the results of the first approxi- 
mation by estimating the variations in flow conditions in the 
0-direction between blade elements on the surface of revolution. 

Assumptions . - The major assumption of the second approximation is 
that, for high-solidity blade-element profiles with minimum radii of 
curvature at least approximately equal to the channel width between 
profiles (except at the nose and tail) , flow conditions and products of 
flow conditions deviate only moderately (and smoothly) from a linear 
variation in the 0-direction across the channel between blade elements. 
Therefore, integrals involving the variation in the flow conditions 
(Q and p, far example) in the 0-direction can be evaluated by numeri- 
cal methods using two-strip integration formulas (that is, ass uming 
parabolic variations in the integrand) . 

In reference 13 it is shown that, for a high-solidity two- 
dimensional cascade with large turning angle and with a minimum radius 
of curvature along the blade-element profile (exclusive of the nose 
and tail) somewhat less than the channel width, the variations in Q, 

Q cos p, and Q sin P can be represented in the 0-direction by a 
parabola. Also, for curved channels constructed from concentric stream- 
lines of a potential vortex it Is easily shown that, if the smaller 
radius of the channel walls is equal to the channel width, the velocity 
at midchannel is only U percent less than that given by a linear 
variation across the channel. Thus far this type of channel, if the 
minimum radius of curvature is at least approximately equal to the 
channel width, the velocity distribution deviates only moderately (and 
smoothly) from a linear variation and may therefore be represented by 
a parabola. 

Flow field . - In the second approximation, the flow field on the 
surface of revolution is limi ted, as before, by the n- and t-boundaries 
(fig. 8) . The flow direction P and the velocity Q vary in the 
©-direction across the channel between blade elements, so that for 
constant M 


P = p(X) 
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and 


Q = Q(X) 

where X indicates position across the channel in the ©-direction in 
percent of total distance between surfaces -of adjacent blade elements 
(fig. 8) 


©l-© 0 

Each streamline in the channel between blade elements is designated by 
the value of the stream function ^ along the streamline where, for 
convenience, the stream function is defined as the ratio of flow rate 
between the streamline and the left side of the channel, when viewed 
in the direction of %£, to total flow rate between blade elements. 
Thus, the 0-streamline, as indicated by the zero value of ljr, flows 
along the left channel wall (fig. 8); the 1.0-streamline flows along 
the right channel wall; and the 0.5-streaml i ne equally divides the flow 
between blade elements, numerical subscripts on p, X, and the 
velocity Q refer to the streamlines at which these quantities are 
considered. Thus, Po.5 ^0.5 refer to the values of p and X 

along the 0.5 -streamline (fig. 8). 

Outline of method . - An inspection of figure 8 shows that the 
blade-element profile on the specified surface of revolution (with 
specified variation in H*) is completely described by the solidity a 
and by 



x o.5 = Xo.sC 11 ) 

(31) 


Po.5 = ^0.5^ 

(32) 

and 

5 = S(M) 

(33) 

These quantities 

(d, Xq.5, Po. 5> 311(1 ' 8 ) 0311 8e 

determined (approxi- 

mately) from the 

prescribed velocities Qq and 

and from the pre- 


scribed upstream and downstream conditions in a manner outlined as 
follows: 
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(a) The solidity a is the same as determined in the first 
approximation . 

(b) The flow directions (approximate) on the blade surface p^ 
and pJ^ are obt ain ed at arbitrarily specified values of M from the 

blade-element profile deter min ed by the first approximation or by the 
previous cycle in the iteration calculations to be used in the second 
approximation. The superscript prime ( 1 ) indicates values obtained 
from the first approximation or freon, the previous cycle in the itera- 
tion calculations. 

(c) The ratio Xq .5 is approximately determined at specified 
values of M from continuity considerations using the prescribed values 
of Qo and Qp and the values of p^ and p^ obtained by (b) . 

(d) The velocity 0o,5 at X^g is estimated at specified values 
of M from considerations of absolute irrotational motion. 

(e) The angle Po.5 i s approximately determined at specified 

values of M from the known value of (Qs)av (equation (22)) assuming 
(see section Assumptions) that integrands consisting of the product of 
various flow conditions can be represented by parabolas that pass 
through known (or estimated) values of Q, P’, and p/p Q at X 
equal to 0, Xq. 5 , and 1.0. 

(f) The channel-width ratio 8 is approximately determined at 
specified values of M from continuity considerations assuming again 
that integrands consisting of the product of various flow conditions 
can be represented by a parabola. 

(g) From the values obtained for X^g, Po.5> a 

element profile is determined so that new values for P^ and p.£ can 
be obtained in (a) and the process repeated until the values of Xq. 5 > 
pQ.g., and 8 converge. 

Flow directions Pq and P£. - In order to obtain the flow direc- 
tion Pq along the blade - element profile determined by the first 

approximation (or by the previous cycle in the iteration calculations 
of the second approx im ation) , equation (30) is written 


e 6 


= e ! 

0.5 


- X 0.5 


( 0 i-e«) 
v 1 o 



KACA TN 2408 


19 


or, from equations (ll) 


and ( 12 ) , 


6 b 


~ 0 O.5 ” X 0.5 


8 « 

Rrpd 


which, after differentiating with respect to M and multiplying "by R, 
"becomes 


ton ?■ = tan f ^ 5 


R* afoo.s S') 
q dM 


(34) 


In like mann er. 


tan ^ = tan P * <5 



(35) 


Equations (34) and (35) determine 0 q and 0|. For the first cycle in 
the second approximation, ^ 6.5 i s equal to 0.5 (equivalent to the 
assumption in the first approximation that the shape of the mean (0.5) 
streamline is the same as the "blade-element camber li n e), and 5* and 
0 ^ 5 are equal to 8 and 0 from the first approximation. Far suc- 
ceeding cyeles in the iteration calculations of the second approxima- 
tion, 0 ^ 5 X^ 5 , and 5’ are obtained from the preceding cycle. 

Ratio X 0 . 5 . - From the definition of the stream function ijr. 


t 


r -£- Q cos 0 dX 
P o 

~ r i 

Q cos 0 dX 

Jo P o 


(36) 


For the 0 . 5 -streamline, i|r equals 0.5 and X 
with 0 equal to 0 1 equation (36) becomes 


equals Xq. 5 > so that 
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0.5 


pxo.s 


Q cos p* dZ 


pl.O 

^0 


■£- Q cos p* dZ 
p o 


(36a) 


If the integrand (t£- Q cos p*) in equation (36a) deviates only 

P ° / 
moderately from a linear variation in the Z-direction (see section 

Assumptions ) , it can he shown (appendix B) that, if a linear varia- 
tion in' the integrand is assumed in equation (36a), the value of Zq .5 
obt ain ed from equation (36a) is approximately equ al to the value 
obtained by representing the integrand by a parabola. That is, for 
moderate deviations of the integrand from a linear variation, the value 
of Zq 5 depends more on the magnitude of the integrand at Z 
equals 0 and 1.0 than on the magnitude of the deviation. If a 

linear variation in (■£- Q cos P*) with Z is assumed, equation (36a) 

p o 

becomes 


0.5 


k Xq.5 + \ Xq.s 2 
k + i 


or 


Xo 


.5 = - k ±yk 2 + k + I 


(37) 


where 


k = 


fe Q °° 5P ')n 

(4 Q C0B “'X - (£ Q 008 P 'Jo 


(37a) 
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The value of Xq.5 value of M. is determined by equa- 

tion (37). These values of Xq .5 can be used to compute new values 
of Pq and by equations (34) and (35), and these new values 

of {J^ and result in new values of Iq ^ from equation (37). This 

iterative process converges rapidly and, in general, the second value 
of Xq 5 obt ain ed from equation (37) is sufficiently accurate. 

Velocity Q0.5. - The velocity Qo.5 is given by the following 
equation (obtained from equation (C8) in appendix C): 


«0.5“ [% + ! ( E + f) + f] ^ < AX 0.5> -!( E+ ! + FX 0.5) ’I 

(38) 

where A is obtained by trial and error from 


where 


% - Qq exp (A) 

+ C 

1 - 

exp (A) 1 

A J 

( e + 1) C 1 

exp 

A 


+ x 


(C9) 


' B * 008 P 0.5 Bln P 6.5 

(C6) 

6* _* TT „ t sin a 

a R U T cos P Q>5 ^ 

(C7) 

E ^ 
E dM 

(C3a) 

<% iQo 

F “ dM " dM 

(C3b) 
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For small values of A, the numerical solutions of equations (38) 
and (C9) involve small differences of large numbers. In order to eli- 
minate this condition, the equations can he expanded in the following 
series form: 

Qo >5 = Qq exp (A X q>5 ) + BE Z Q ^ 5 + C X Q>5 + 

_.s( X 0.5 2 AX 0.5 5 a2 ^.s 4 \ 

(BEA + BF + CA) ( 2 , + 3 , + 4 , + . . .1 

and 

_ [% - Op eg (A)] - C - CA ( Jj. + A. + + . .) 

E + (a + D + + • • •) 

in which forms the calculations are mere easily carried out for small 
values of A. 

Flov direction p Q j5 . - At any position (X,M) in the channel 
between blade elements, the flow direction 8 (determined by the second 
approximation) is related to the uncorrected flow direction P’ (deter- 
mined by the blade-element profile of the first approximation, or by the 
profile of the previous cycle in the iteration calculations of the 
second approximation) by 


(38a) 


(C9a) 


P = p* + AP (39) 

where Ap is the -correction added to p 1 . In general, for a given 
value of M, the correction Ap varies with 0, but it can be shown 
(appendix D) that in the direction of 0, if (tan P - tan P 0 . 5 ) is 
small (a condition that is approached in high-solidity blades with 
gradual variations in blade-element-profile thickness) and if APq .5 is 
small, then 


Ap ~ ^Pq.5 

The correction Ap is therefore assumed constant across the channel 
between blade, elements. Thus, for example, equation (39) becomes 
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P 0.5 = P0.5 +A P 


(40) 


The correction Ap is determined from considerations of the 
average relative tangential velocity, -which is defined "by 


Ol 


(^fPav - 


0 P o 


7 ^- Q 2 sin P cos p dX 


AL 


■£- Q cos p dX 


From equation (39) -with Ap assumed to he s mall , 

cos p « cos P* - Ap sin p 1 
sin p » sin p* + Ap cos p* 

and 


sin P cos p - — sin 2P * + Ap cos 2p * 

u 


From equations (41) and (42), 


(%)avj^ ^ Q cos P' dX - ^ Q 2 sin 2p« « 

pi pi 

Q 2 cos 2P ' dX + (0e) aT #- Q sin p* dX 

JO ° Jo ° 


(41) 


(42a) 

(42b) 

(42c) 


(43) 


Let I he any one of the integrands in equation (43) ; then if I 
varies in a parabolic manner (see section Assumptions ) across the 
channel between blade elements. 


= % + ®b x + e c 


I 0 = e a 

*0.5 = «% + ®b Xq. 5 + e c ^O.S 2 


so that 



24 


NACA 031 2408 


and 


h = e a + % + e c 


Therefore, 

I ( Z *0-5 - 1 \ ( 1 \ /3 Xp.5 - 2 \ 

_°V SXo.5 )’ °- 5 (,3Xo. 5 (Xo. 5 -l)J ^3 20.5 - 3; 

. (44) 

(With Xq .5 equal to 0.5, equation (44) reduces to Simpson’s one-third 
rule for numerical integration.) 



The direction Po.5 4 s therefore given by equation (40) in which 
AS is obtained from equation (43) where the integrals are determined hy 
equation (44) . (The values of Q, p/ p Q , and p ’ contained in Iq, 

Iq. 5 ^ and I]_ of equation (44) are known. Far the first cycle of the 
second approximation, p^ is equal to p obtained hy the first 
approximation and far succeeding cycles of the second approximation, 

Pq 5 is obtained from the preceding cycle. ) 


Channel -width ratio S . - From continuity considerations of the 
fluid flow across the strip between channel walls in figure 7, 


f 0 l 

~ = [cH* (AH)J (cR* Bj P 0 a 0 #- Q cos p d£? 

L -J - o 


■which, from equations (l2), (19), (20), and (30), becomes 


cp = 6 H*R* 


f Q cos p dX 
o 


or, from equations (39) and (42a), 
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8 = 



Q cos p ‘ 
K o 


CP 


dX - Ap Q sin p* ax 

Jo ° 


C^5) 


Equation (45) determines tlie channel -wi dth ratio. The integrals in 
equation (45) are given by equation (44) . 


The values of Po.5 sad- $ obtained by the first cycle of the 
second approximation can be used to determine new values of p^, Pj^, 

Xq 5 , and Q 0.5 from which new and presumably better values of Pq .5 
and 8 can be determined. This iterative process can be continued 
until the values of Po.5 8 converge. In general, the iterative 

process converges rapidly and, in many cases, only the first cycle of 
the second approximation is necessary. 


Blade- element-profile coordinates . - The value of O has been 
determined, and R*, H*, X 0 . 5 .J Po.5> and. ^ 8X6 know 11 functions 

of M so that the design of the blade element is complete. The coordi- 
nates of the blade-element profile on the surface of revolution are 
obtained as indicated in figure 8 by (l) computing the shape of the 
0.5-streamline, which is obtained from Po.5 = Pq. 5 ^ > ( 2 ) distri- 

buting the channel width on either side of the 0. 5-streamline according 
to Xq .5 = Xq. 5 (m), and ( 3 ) obtaining the blade profile from the 
channel shape and from the spacing of the channels, which spacing is 
determined by cr. 

The shape of the 0.5-streamline is given by R and 9 as func- 
tions of M. Because R (or R*) is a known (specified) function 
of M, only (R0 )q. 5 remains to be determined as a function of M. 
From equation (28) , which assumes R0 is zero at M = 0, 


(»>0.5 



(tan Pq 5 + ® sin a ) <3M 


(46) 


Equation (46), together with the specified variation in R, determines 
the shape of the 0.5-stre amlin e on the surface of revolution. 
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The channel width is — R* of which Xa c -H* is distributed 

a a 

e 

to the left of the 0.5-strea ml iae and (l - X 0 . 5 ) — R* is distributed 

to the right {when viewed in the direction of Qy, fig. 8 ). The 
expression far the channel width is dimensionless and is expressed, as 
usual, in units of characteristic length c. 

The blade-element profile results when the adjacent channel is 
placed with the spacing between 0.5-streamlines equal to R*/cf 
{fig. 8 ) . The blade-element profile is closed by rounding-off the noBe 
and tail in a mann er guided by experience. 


APPLICATION CEF DESIGN METHODS 

The application of the first and second approximations of .the 
blade-element-design method is considered in this section. 

Types of blading . - The blade-element -design method developed in 
this report applies to compressible or incompressible flow in stators or 
rotors of radial-, axial-, or mixed-flow c empress or s, turbines, or two- 
dimensional cascades. The design of any particular blade element 
depends on which of the preceding categories describe the application 
of the blade element. Each category determines the magnitude or sign 
of various parameters that appear in the equations developed for the 
design method. 

The blade element is first classified according to the compres- 
sibility of the fluid for which it was designed. For compressible 
fluids, the density ratio p/p Q varies according to equation (lO)j 

for incompressible fluids, the density ratio is constant and equal to 1 . 

The blade element is next classified as a stater or rotor. For 
stators the angular velocity cd Is zero so that the tip speed Uj is 
zero (equation (4)); for rotors the tip speed is positive. (Note that 
the angular velocity is always positive (right-hand rule) about the 
Z-axis (fig. 4) . ) 

The blade element is next classified as inflow, axial, flow, or 
outflow. For inflow blades (fig. 4(a)), the radius R decreases in 
the direction of M; thus from equation (la) the angle a is less than 
zero. For axial -flow blades, a approaches zero. For outflow blades 
(fig. 4(b)), the radius R increases in the direction of M so that a 
Is greater than zero. 
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The "blade element; is also classified as a compressor, turbine, or 
two-dimensional-cascade "blade. In the case of rotor blades, Qi is 
greater than Qq for compressors, and Qq_ is less than Qq for tur- 
bines. For plane two-dimensional cascade blades, the surface of 
revolution becomes a flat plane, and 


H* = 1.0 
R* = 1.0 


R =oo 
dB = 0 


Rde = an > 


(47) 


a = 0 


u T = 0 


J 


where N is the coordinate distance measured normal to M on the flat 
plane of the two-dimensional cascade. 

Prescribed upstream and downstream flow conditions . - The designer 
of blade elements is generally given the upstream and downstream flow 
conditions for each blade element. These flow conditions are: 

(l) specified surface of revolution R* = R*(M), (2) specified height 

ratio H* = H*(m) , (3) tip stream and downstream whirl ratios 

and Aj), and (4) the flow coefficient Cp (equation (l9) ) . 

For axial-flow blade elements, R* and H* are considered con- 
stant upstream and downstream of the blade-element row so that the 
quantities Qq, Pq, Qq, and Pq are constant. It is therefore cus- 
tomary to specify any three of these quantities (instead of Aq, Aq, 
and cp); the fourth quantity is then determined from continuity con- 
siderations (equation (26)). For radial- or mixed-flow blade elements, 
however, R* and H* generally vary upstream and downstream of the 
blade-element row so that (Q and P)-q and (Q and p) D vary and it 
is more convenient to specify the upstream and downstream conditions by 
the quantities Aq, Aq, and cp, which are always constant. 

Prescribed velocities along blade-element profile . - In this 
blade-element-design method, the prescribed velocities Qq and Qq on 
the blade-element surfaces are specified as functi ons of the meridional 
coordinate M, rather than as functions of the arc length S, along 
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the blade -element profile. Any relation "between Qq ahh Qp and M 
can "be prescribed, "but, if a physical l y acceptable blade-element profile 
is desired, and if gradual variations in blade-element ^-profile thickness 
together with a minimum radius of curvature at least approximately equal 
to the channel width is desired in order to satisfy the assumptions of 
the design method, the prescribed relation must exhibit the following 
features: 

(l) In order to avoid negative blade thicknesses, which result 
when the channel width in the circumferential direction is greater than 
the blade spacing (that is, when 5 >1.0), the average velocity Qgy, 
defined by 


Qav = | foo + Ql) 


(24) 


must not be too low. As indicated by equation (25) , the minimum allow- 
able value of Qgy depends, among other things, on the value of p, 
which unfortunately is not known until the blade element has been 
designed. 


(2) In order to obtain desirable blade-element thicknesses 



at the n- and t-boundaries (fig. 8, for example), the specified values 
of (Qa V ) n and (Qav^t must be approximately equal to the values given 
by equation (25) for specified values of 5 and estimated (not yet 
known) values of p. These values of Qgy must be adjusted later if 


they result in undesirable blade-element thicknesses in the final cal- 
culations. 


(3) In order to 
blade thickness near the nose. 


obtain a physically desirable rate of increase of 
b-Qav'' 


(W* v) 

\ <3M / r 


at M = 0 must have a suitable 


dM /ji 

positive value, which is best determined from experience and which 
depends on the unknown solidity (to be determined), on the unknown 
variation in P with M (to be determined) , and on the prescribed 
variation in a H with M. Lik ewise, to obtain a desirable rate of 

/dQav\ 

V <3M J ± 


decrease of blade thickness near the tail. 


at M = 1.0 must 


have a suitable low, perhaps negative, value, which also depends on 
the solidity and on the variation in p and R*H* with M. 



NACA Iff 2408 


29 


(4) In order to obtain the fewest "blades (that is, lowest solidity) 
possible, consistent with other requirements and limitations, it is 
desirable that the blade loading AQ, defined by 

AQ = Q-l - Qq (48) 

be as large as possible over as great a range of M as possible. This 
fact is indicated by equation (21) , because AP increases with AQ. 
Thus, it is generally desirable to prescribe blade loading even at the 
nose and tall (n- and t-boundaries, fig. 8). The blade is then assumed 
to load and unload over the regions that correspond to rounding-off of 
the nose and tail, respectively. The loading (AQ) n and (AQ)-t must 
not, however, be too large; otherwise assumptions of the design method 
are violated. For example, if AQ is large in the vicinity of the 
nose or tail, this sudden loading or unloading of the blade can be 

expected to cause p, Q, and — to deviate greatly (and in a sudden 

dM 

manner) from a linear variation in the 0-direction (which deviation 
violates the assumption of the blade-element design method that the 
deviation is smooth and moderate) . The maxi-mum loading (AQ) m should 
occur at some value of M between, but not close to, the nose and 
tail (0 < M < 1.0). 

(5) The prescribed variation in Qo and Qj_ with M should be 

smooth with gradual variations in ^2. and Sudden variations in 

dM dM 

dQo dQi 

and result in small radii of curvature locally on the blade - 

dM dM 

element profile (such as occur at the nose and tail, for example) and 

cause p, Q, and — to deviate greatly (and in a sudden mariner) from 

dM 

a linear variation in the 0 -directi on. 

One possible way of incorporating the preceding considerations in 
the prescribed distributions of Qq and Q^ is given by 


<fl w 

i 

II 

(49a) 

% = Q av + f 

(49b) 
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■where 

Qav = ^ ^ M + d e M 2 + M 3 (50a) 

and 

AQ = d a + h M + j c m2 + % M3 (50b) 

The coefficients d and j in equations (50a) and (50b) are obtained 
from the following conditions: 


so that 


= teaA at M = 0 

Qav = C^av^t at M = 1.0 


dQav 

= (Wav) 

at 

M 

= 0 

dM 

V CM A 





. ( dC M 

at 

M 

= 1 

dM 





AQ = 

= N n 

at 

M 

= 0 

AQ = 

= Nt 

at 

M 

= 1 

AQ = 

= (AQ) m 

at 

M 

= M m 

d(AQ) 

dM 

* 0 

at 

M 

= % 
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and 


^a ~ COav^n 

, _ (^a.v\ 

^ \ < 3 M y n 

^ = 3 [(Qav)t - CQav)n] “ (^) t " 2 (^f) n 
«» - (^)n + (^)t - 2 [Wt - COn] 


(51) 


J a = (AQ) n 


3b 


1% 2 (l-Mj 2 


" C^)n (Mm 4 - 3V 2 + 2MJ - 


<3a “ 


(AQ) m (S^ 2 - 2Mj + (AQ) t (m/) 


(AQ) n (ZmJ 5 - a^ 2 + 1) + 


% 2 ( 1 -Mm ) 2 


3d 


(AQ) m (3^ 2 - 1) - (AQ) t (2^3) 


“ (^)n (Mm 2 - 2^ + !) - 


Mm. 2 


(AQ) m (2^-1) + (AQ) t (l^ 2 ) 


J 


(52) 


Equations (49a) and (49b) are a suggested means of prescribing 
satisfactory distributions of Qq and Q^; any other method that 
results in acceptable blade-element profiles can be used. 
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Outline of numerical procedure . - The following outline of the 
numerical procedure includes both the first and second approximations. 
For convenience, the outline is divided into four parts: (i) Speci- 

fied Conditions, (il) Preliminary Calculations, (HI) First Approxima- 
tion, and (iV) Second Approximation: 

I - Specified Conditions 

(l) Specify R = E(z), which determines the specified surface of 
revolution and from which is obtained the constant Erp and the rela- 
tions 


R* = R*(M) 


and 


a = a(M) 


(2) Specify H* = H*(M) . 

(3) Specify Ag, Ag, and cp, which determine the flow conditions 

upstream and downstream of the blade elements. (For axial-flow 
blade elements, it is customary to specify (instead) any three of the 
quantities Qg, Pg, Op, and Pd; the fourth quantity is then deter- 
mined by equation (26) in which is given by equation (10) . The 

quantities Ag, Ag, and cp are then determined by equations (8) 
and (26)). 

(4) Specify Uqi, which is zero for stators. 

( 5 ) Specify Qg = Qg(M) and Qg = Qg(M) . 

(Use equations (49a) and (49b) or any other arbitrary method. These 
prescribed values of Qg aa<I ®1 may have to be adjusted later to 
result in physically acceptable blade-element profiles on the surface 
of revolution.) 

II - Preliminary Calculations 

(1) Compute AP as a function of M from equations (9) and (16). 

(2) Compute the solidity a from equation (2l) by numerical inte- 
gration. 
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(3) Compute number of "blades b from 


. 2it 

b ‘M 


where A 9 is obtained from equation (ll) . If the number of blades is 
not an integer, the solidity a must be changed by adjusting any or 

all of cp, Aj), Ay, and f H*B*(AP) <3M, as indicated by equa- 

J 0 

tion (21) . 

(4) Compute (Qe)av as a function of M from equation (22) by 
numerical integration. 


Ill - First Approximation 

(1) Compute Qqv as a function of M from equation (24) . 

(2) Compute (3 from equation (23) . 

(3) Compute 8 from equation (25) . 

(if the coordinates of the blade-element profile on the surface of 
revolution are to be determined after the first approximation (that is, 
if the second approximation is not used), continue with the remaining 
steps in this section; otherwise proceed immediately from here to the 
next section.) 

J 

(4) The R,6 coordinates of the blade camber line on the surface 
of revolution are given by the specified relation between R (or R*) 
and M, and by R0 as a function of M, which is obtained from equa- 
tion (28) by numerical integration. 

(5) The blade thickness in the circumferential direction is equally 
distributed on either side of the blade camber line and is given by 
equation (29). 

(6) The angular blade spacing A 9 is obtained from the solidity, 
according to equation (ll) . 

(7) The nose and tail of the blade are rounded off (fig. 6) in a 
mann er guided by experience. 
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( 8 ) For plane two-dimensional cascades tlie M,0 coordinates are 
replaced by M,IT coordinates where dN replaced the product Rd0. 
The solidity cr, defined by equation (ll) , thus becomes equal to the 
reciprocal of the blade spacing (which spacing is expressed in units 
of the characteristic length c, so that the solidity is equal to 
the chord divided by the blade spacing) . 

XV - Second Approximation 

The second approximation (part IV) follows after step (3) of the 
first approximation (part in) . 

(1) The angles 8 q and 8 ^ are obtained from, equations (34) 

and (35) with Xq 5 assumed to equal 0.5, and with 81 and 8 *’ 

* u .0 

equal to 8 and 8 , respectively, from the first approximation. 

(2) The ratio Xq^ is obtained from equation (37). 

(3) If desired, the angles 8 q and 8 ^ can be recomputed from 
equations (34) and (35) using Xq .5 from step (2). This iterative 
process can be continued, although the process generally converges 
rapidly enough so that only one step in addition to step ( 2 ) is neces- 
sary. 


(4) Compute the corrected flow direction 8 q 5 from equations (39) 
and (43) . The integrals in equation (43) are evaluated by equation (44) 
in which the velocity Qo ,5 is given by equation (38) or (38a). To 
solve these equations initially, the valueB of 8 q 5 and S’ are the 

8 and 8 obtained by the first approximation. 

( 5 ) Compute the corrected channel-width ratio 8 by equa- 
tion (45) . 

(6) The corrected values of Pq 5 ^ are obtained initially 

from quantities that are based upon 8 and 8 determined by the first 
approximation. If these quantities (8qj 8-[> Xq and Qq g, for 

example) are recomputed using the corrected values of Po.5 a^i 8, 
new, and presumably better, values of 8 q g and 8 are obtained. 

This iterative process can be continued until the values of P 0.5 

and 8 converge. In general, the process converges rapidly and, in 
many cases, only the first cycle of the second approximation is neces- 
sary. 
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( 7 ) The shape of the 0. 5-streamline on the surface of revolution 

is given by the specified relation between B (or B*) and M, and by 

(B0)o. 5 = f ( M ) 

which is obtained from equation (46) by numerical integration. 

(8) The channel width — E* is distributed on either side of the 

o 

0.5 -streamline according to the ratio Xq^, which is obtained from 
equation (37) using the final corrected values of 3o.5 and 5. 

(9) The blade-element profile results when the adjacent channel is 

placed with the spacing between 0.5-streamlines equal to — (fig. 8). 

cr 

(10) The blade-element profile is closed by rounding off the nose 
and tail in a manner guided by experience. 

Boundary-layer corrections . - The blade profiles obtained by the 
first and second approximations are based upon the assumption of a non- 
viscous fluid. For real, viscous fluids a boundary layer, within which 
the viscosity cannot be neglected, builds up along the blade surfaces 
and displaces the potential flow. If the prescribed velocities for 
which the blade was designed do not have large enough decelerations to 
cause separation of the boundary layer, the thickness of the boundary 
layer can be corrected for, at least partially, by machining from the 
blade profile an amount equal to the displacement thickness of the 
boundary layer, which thickness can be estimated from boundary-layer 
theory. 


NUMEBICAL EXAMPLES 

Two numerical examples are presented. The first example is for 
compressible flow between blade elements of a plane two-dimensional cas- 
cade with prescribed velocities along the blade-element profile. The 
second example is for compressible flow between blade elements of a 
mix ed-flow impeller with prescribed velocities along the blade-element 
profile. In both examp les, the design method has been checked by com- 
paring the prescribed velocity distribution with the velocity distri- 
bution obtained by stream-filament methods on the resulting blade- 
element design. 
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First numerical Example 

The first numerical example is the design of a "blade -element pro- 
file for a plane two-dimensional cascade that accelerates the fluid by 
turning it from 6 equals zero upstream of the cascade to (5 equals 
60° downstream. 

Prescribed conditions . - For any plane two-dimensional cascade the 
conditions given hy equation (47) apply. Also, for the first numerical 
example the following flow conditions upstream and downstream of the 
cascade are given: 


Qd = 0.750 
= 60° 
Pg= 0° 


so that, from equation (26), 


Qq = 0.290 

and 

Cp = 0.278 

Prescribed velocity distribution . - The prescribed velocity dis- 
tribution was determined by the following quantities, the magnitudes 
of which were based on considerations previously outlined in this 
report: 


(Qav^n — 0-348 
(Qav)t = 0-751 


by trial and error to obtain 
satisfactory values for 6 n and 6^ 



0.444 




chosen so that 



0 
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(AQ) n = 0.200 
(AQ) t = 0.100 
(AQ) m = 0.400 
= 0.450 

so that, from equations (50a), (50b), (5l), and (52), 

Qg^ = 0.348 + 0.444 M + 0.371 M 2 - 0.412 M 3 

AQ = 0.200 + 0.888 M - 0.984 M 2 - 0.004 M 3 

from which Qq and Qq_ are obtained by equations (49a) and (49b) . The 
velocity distribution is given in table I and is platted in figure 9. 
(Sonic velocity occurs at Q equal to 0.915.) As a result of this 
velocity distribution the solidity a is equal to 1.215 (from equa- 
tion (21)). The pressure difference AP and the average relative 
tangential velocity (Qe) a v tabulated as functions of M in 

table I. 

Results . - The values of 0 and 5 that determine the blade- 
element profile for the first approximation are given in table I 
together with the values of Pq g, 8, and Xg^g that determine the 
blade profile for the second approximation. The blade-element profile 
obtained by the second approximation is plotted in figure 10 together 
with the streamlines (obtained from equation (36) ) and velocity poten- 
tial lines used by the stream-filament method to estimate velocities 
along the blade-element profile. (The blade nose and tail are rounded 
off in an arbitrary manner.) The velocities obtained by the stream- 
filament method are compared in figure 9 with the prescribed veloci- 
ties. The comparison is considered satisfactory, although the design 
method is not recommended for lower solidities. 

The resulting blade-element profiles for the first and second 
approximations are compared in figure 11. The two profiles are sig- 
nificantly different and it is concluded that for blade-element soli- 
dities as low as that for the first numerical example (cr = 1.215), the 
second approximation is required. 



38 


MCA TN 2408 


Second Numerical Example 

The second numerical example is the design of a "blade-element 
profile for the Impeller of a mixed-flow centrifugal compressor with 
prewhirl (Ay) and with forward-curved blades at the im peller tip. 

Surface of revolution . - Far the mixed-flow impeller of the second 
numerical example, the surface of revolution is generated by rotating 
the circular arc shown in figure 12 about the axis of the impeller 
(Z-axis) . At the blade-element nose (n-boundary) 


On = 0 


= 0.5 Rrji 


and at the blade-element tail (t -boundary) 


oj. = 45° 


R^; — 1.0 Rjl 

so that the circular arc subtends an angle of 45° and has a radius Y 
equal to 1.707 Rrp. The length of the circular arc, or blade-element 
chord, is unity, so that (from fig. 12) 

Y y = 1.0 
4 

from which R^ is equal to 0.746. 

Transformed coordinates . - For convenience in the stream-filament 
check the following transformation of coordinates, which conformally 
transforms the surface of revolution into a flat plane, has been 
introduced: 


and 




(54a) 


(54b) 



HACA TIT 2408 


39 


where 


dM 

cLL = = 


dR 


R R sin a 


and . 


AL = 
u 


dL = 

Ln ^ 


pRfc 


dR 


R sin a 


*n 


For the surface of revolution in the second numerical example. 


(55a) 


(55b) 


1 _ Y 

sin a I p ' p 

V Y - (Y + R^ - R) 2 

so that equation (55b) integrates to give 


(56) 


r 


AL = 


^ (21 + Rj 


Jt , . -1 

— + sin 

u 


Rt (Y + Ra) ~ Rq (ZY + R]^) 


>= 2.096 


Also, equation (54a) , combined with equations (55a) and (56) , inte- 
grates to give 


AR 


dR 


R sin a 


L* = 




AL 


$ + sin -1 
2 

_ R (Y + Rq) - Rjj (2Y + Rjj) 

L RY 

— + sin“ Y 
2 

Rfc (Y + - Rn (ZY + Rj" 

L RuY 


(57) 


Equation (57) deter min es R (or R*) as a function of L*. The quan- 
tity L* varies from 0 at the n-boundary (fig. 13(a)) to 1.0 at the 
t-boundary, as does the coordinate M. Equations (56) and (57) deter- 
mine a as a function of L*. The quantities R* and a are tabu- 
lated as functions of L* in table II. 
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Lines of constant L * and 0* on the surface of revolution and 
on the flat transformed plane are shown in figures 13(a) and 13(h), 
respectively. The equations in this report were developed for surfaces 
of revolution (<3M and d0 coordinates), hut, if dM is replaced 
hy b(al) <3L , the equations also apply to the flat transformed 
surfaces . 

Prescribed conditions . - In addition to the surface of revolution 
which has a l ready "been described, the following conditions are pre- 
scribed: 

H* = 15-14R* 

Urp = 1.5 

= 0.2 
= 2.0 

and 

cp = 1.05 

These prescribed conditions were selected to result in (approximately): 

(Qa V )n = 0-50 
Pn = - 45° 

Sn = 0.85 
(Qav^t = 0.90 
5 t = 0.95 

Prescribed velocity distribution . - The prescribed velocity dis- 
tribution for the second numerical example was obtained from equa- 
tions (50a) and (50b) with x instead of M as the independent vari- 
able, where 


x = 2R* - 1 


(58) 
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The prescribed conditions an the velocity distribution are: 

(Qav)n =0-5 

(Qav)t =0.9 

chosen so that 

(ttX " °- 55 } oiOBen 30 that (fr)t = 0 

(AQ) n = 0.1 
(AQ) t =0.2 
(AQ) m = 0.5 
% = 0.5 

so that, from equations (50a), (50b), (5l), and (52), 

Qav = 0.50 + 0.35x - 0.05x 2 + 0.10x3 
AQ = 0.10 + 1.70x - 2.00x2 + 0.40x 3 
from which Qq and Qq_ are obtained by equations (49a) and (49b) . 

The resulting velocity distribution as a function of L* is given 
in table II and is plotted in figure 14. As a result of this velocity 
distribution the solidity a is equal to 2.561 (from equation (21)), 
and the number of blades b is 12 (from equation (53) ) . In order to 
obtain this even number of blades, the downstream whirl Ap was 
adjusted (as indicated by equation (2l)) from the prescribed value 
of 2.000 to 1.939. The pressure difference AP and the average rela- 
tive tangential velocity (Qg) a v 31,0 tabulated as functions of L* 
in table II. 
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Results . - The values of £J, 5, and R* that determine the 

"blade-element profile for the first approximation are given as func- 
tions of L* in table II together with the values of Po.5> 8, Xq, 5 > 

and R* that determine the blade-element profile for the second 
approximation. The blade-element profile obtained by the second 
approximation is plotted on the surface of revolution in figure 13(a) 
and on the flat transformed L*S*-plane in figure 13(b). (The nose 
and tail of the blade-element are rounded off in an arbitrary manner.) 
The blades are curved forward in the direction of rotation at the 
impeller tip, because the prescribed value of Ap is greater than Urp. 

The variation in thickness of the blade-element profile at values 
of L* toward the impeller tip is considered rapid because, for the 
relatively large values of {3 that exist at these values of L*, the 
difference (tan - tan Pq) is large, which condition makes inac- 
curate the assumption that Af3 is constant circumferentially across 
the channel between blade elements (see section on flow' direc- 
tion Po.5)* The rapid variation in thickness can and should be eli- 
minated by different prescribed variations in H* and Qq^ or both 
with L*. 

The streamlines and their normals, used by the stream-filament 
method to estimate velocities along the blade-element surfaces, are 
plotted on the L*0* -plane in figure 13(b) . 

The velocities obtained by the stream-filament method on the 
L*0* -plane are compared in figure 14 with the prescribed velocities 
as a function of L*. The comparison is excellent, except at large 
values of L* where the rapid variation in blade-element-profile 
thickness and the small radii of curvature along the blade-element 
profile at the tip (fig. 13(a)) violate assumptions of the design 
method. Thus the design method is limited to prescribed velocities 
that result in blade-element profiles with gradual variations in 
thickness and with minimum radii of curvature at least approximately 
equal to the channel width between profiles. 

The resulting blade-element profiles for the first and second 
approximations are compared in figure 15. The two profiles sire quite 
similar, except near the tail (where neither approximation is satis- 
factory because of the rapid variation in blade-element-profile thick- 
ness and the small, radii of curvature at the blade-element tip) . 
Although the difference in profiles is small, the difference may be 
significant because of the small blade spacing. However, at least for 
rough calculations, the first approximation appears to be satisfactory 
for centrifugal Impellers with blade-element solidities as high as 
that of the second numerical example (0 » 2.5). 
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SUMMAEff OF RESULTS AMD CQNCLUSIOUS 

An approximate blade-element- design method is developed for com- 
pressible or incompressible ncnviscous flow in high-solidity stators or 
rotors of axial-, radial-, or mixed-flow compressors, turbines, or two- 
dimensional cascades. The method is based upon channel -type flow 
between blade elements on a specified surface of revolution that is 
concentric with the axis of the compressor or turbine. The blade ele- 
ment is designed far prescribed velocities along the blade -element pro- 
file as a function of the meridional, distance. The blade-element- 
design method is developed in tiro parts, called the first and second 
approximations, either of which may be used to determine the blade- 
element profile. The first approximation, 'which is simple and rapid, 
assumes that circumferentially across the channel between blade elements 
the velocity is constant and the flow direction is equal to that of the 
blade-element camber line. The second approximation, which is more com- 
plex and less rapid, partially corrects for the variation in velocity 
and for the change in flow direction across the channel between blade 
elements. 

Two numerical examples are presented. The first example is the 
design of a blade-element profile for a plane two-dimensional cascade 
in compressible flow with prescribed velocities along the profile. The 
second example ife the design of a blade element for the impeller of a 
mixed-flow centrifugal compressor. As a result of these numerical 
examples, the following conclusions have been drawn: 

1. The blade-element-design method is limited to prescribed velo- 
cities that result in blade-element profiles with gradual variations 
in thickness and with minimum radii of curvature at least approxi- 
mately equal to the channel width between profiles. (The nose and tail 
regions of the blade-element profile are not considered by this blade- 
element-design method.) 

2. For centrifugal impellers with blade-element solidities as 
high as that of the second numerical example (2.561), the first approxi- 
mation appears to be satisfactory, at least for rough calculations. 

3. The blade-element-design method is not recommended far solidi- 
ties lower than that of the first example (l.215). 

4. Far blade -element solidities as low as that for the first 
example (l.215), the second approximation is required. 


Lewis Flight Propulsion Laboratory, 

Rational Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 26, 1951. 
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APPENDIX A 
SYMBOLS 

The following symbols axe used in this report: 
coefficient, defined By equation (C5) 

characteristic Telocity, absolute stagnation speed of sound 
upstream of Blade -element row where conditions are uniform 
in circumferential direction 

coefficient, defined By equation (C6) 

number of Blade elements (or passages Between Blade elements) 
coefficient, defined By equation (C7) 

characteristic length. Blade chord measured from n-Boundary 
to t -boundary along meridional line on surface of revolution 

coefficient, defined By equation (5l) 

coefficient, defined By equation (C3a) 

coefficient 

exponential, exp (y) = e y 
coefficient, defined By equation (C3b) 

total (Blocked and unblocked) incremental annulus area normal 
to direction of at Blade -element tip, defined By 

equation (20) 

acceleration due to gravity 

incremental height of fluid particle or Blade element, meas- 
ured normal to surface of revolutiop ( fig. 5), (dimension- 
less, expressed in units of characteristic length c) 

height ratio of Blade element or fluid particle on surface 
of revolution (fig. 5), AH/ (ah)^ 

integrand 
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^linear 

I* 

j 

k 

L 

AL 

L* 

M 

N 

P 

AP 

P 

Q 

AQ 

R 

R* 

S 


value of I if Variation in I with X is linear 

integrand ratio, defined "by equation (Bla) 

coefficient, defined by equation (52) 

coefficient , defined by equation (37a) 

transformed coordinate, defined by equation (55a) 

change in L from n-boundary to t -boundary, defined by equa- 
tion (55b) 

(L-Lq) divided by AL (so that L* varies from 0 to 1.0 
along meridional line between n-boundary and t -boundary) , 
defined by equation (54a) 

coordinate distance along meridional line on surface of revo- 
lution (fig. 6, for example) (dimensionless, expressed in 
units of characteristic length e; varies from 0 at n- 
boundary to 1.0 at t-boundary) 

coordinate distance measured normal to ,M on flat plane of 
two-dimensional cascade (dimensionless, expressed in units 
of characteristic length c) 

static -pressure ratio, p/p Q , defined by equation (9) 

change in P circumferentially across channel between blade 
elements, defined by equation (16) 

static (stream) pressure 

relative velocity on surface of revolution (dimensionless, 
expressed in units of characteristic velocity a Q ) (fig. 3) 

change in Q circumferentially across channel between blade 
elements, defined by equation (48) 

cylindrical coordinate (dimensionless, expressed in units of 
characteristic length c) (fig. 3) 

ratio, R/Rj, 

distance along stream lin e on surface of revolution ( dime nsion - 
less, expressed in units of characteristic length c) 

(fig. 3) 
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static (stream) temperature 

tip speed of "blade element, defined "by equation (4) 

total incremental flow rate along surface of revolution with 
incremental annulus height AH 

ratio (circumferential distance from 0 -streamline to point 
in channel divided by channel width in circumferential 
direction), defined by equation (30) 

defined by equation (58) 

radius of meridional ' line on surface of revolution in second 
numerical example (dimensionless, expressed in units of 
characteristic length c) (fig. 12) 

cylindrical coordinate (dimensionless, expressed in units of 
characteristic length c) (fig. 3) 

slope of surface of revolution, equation (l) (fig. 3) 

direction on surface of revolution, equation (2) (fig. 3) 

correction for direction p, defined by equation (39) 

uncorrected flow direction on surface of revolution, defined 
by equation (39) 

ratio of specific heats 

channel-width ratio (circumferential width of channel divided 
by circumferential spacing of blades) , defined by equa- 
tion (12) 

uncorrected channel-width ratio (from first approximation or 
from previous cycle of second approximation) 

cylindrical coordinate (positive about Z-axis according to 
right-hand rule) (fig. 3) 

angular spacing of blades about Z-axis, defined by equa- 
tion (53) 

coordinate 9 divided by AL, defined by equation (54b) 
whirl ratio, defined by equation (8) 
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| integrand ratio, defined "by equation (Bli>) 

p static (stream) density 

cr blade -element solidity, defined by equation (ll) 

cp flow coefficient, defined by equation (19) 

i [r stream function, defined by equation (36) 

a) angular velocity of rotor (in direction of positive 0) 

Subscripts : 

a,b,c,d successive points 

av mass -weighted average value across channel in 0-direction 

between blade elements 

D boundary downstream of blade elements where conditions are 

unif orm in 0-direction 

M component along meridional line on surface of revolution 

m quantities associated with ma xim um value of AQ 

n n-boundary near nose of blade-element profile (fig. 6, for 

example) 

o absolute stagnation condition upstream of blade-element row 

where conditions are uniform in circumferential direction 

R, 0,Z components in positive R-, 0-, and Z-directions, respectively 

T tip (maximum radius) of blade element 

t t -boundary near tail of blade -element profile (fig. 6, for 

example) 

U boundary upstream of blade elements where conditions are 

uniform in circumferential direction 

0 left channel wall, when viewed in direction of Qj/[ (fig* 8) 

0.5 -streamline that equally divides flow between blade 
elements (fig. 8) 


0.5 
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value at X equals 0.5 

right channe l wall^ "when viewed in direction of (fig. 8) 
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APPENDIX B 


VARIATION" H X Q5 
OF INTEGRAND 


FOR VARIOUS PARABOLIC DISTRIBUTIONS 
/^- Q cos p j IN X-DIRECTION 


The ratio Xq ^ is defined "by equation (36a) . This ratio deter- 
mines the position of the mean (0.5-) streamline, -which equally divides 
the flow between two blade elements. The ratio Xq 5 is therefore the 

value of X that eq ually divides the area under the curve of the 


integrand I (equal to — Q cos p 1 ) against X. If I is assumed to 

p o 

vary in a parabolic manner (see assumptions of second approximation) and 
if I has the values Iq at X equals 0, Ii at X equals 0.5, 
and I-, at X equals 1.0, then s' 



where 

* 

I 


i 


in which 

(iiinear^ = 2 ( -^o + ^l^ (Blc) 

IT 

The integrand I in equation (Bl) is dimensionless, having been 
divided by (lunear^i ’ which from equation (Blc) is equal to the value 

s' 

of I at X equals 1/2 if the variation in I with X is linear. 

It was convenient to divide I by (liinear)i ^ ort 3- er to facilitate 

s' 

the comparison of Xq 5 for linear and parabolic variations in I 

with X. The constant I* contained in equation (Bl) is the ratio at 

s' 

X equal to 0.5 of the integrand I x for a parabolic variation in I 

S 


(-^linear^i 


A 

% 


(Bla) 


(Bib) 
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to the integrand 


^linear 


I* 


is a direct measure of 


Ihe 


for a linear variation in I. Thus, 
deviation of the parabolic variation in 


I 2 from a linear variation between the same values of Iq and I 


The 


ratio | (equation (Bib)) is a direct measure of the relative magnitudes 
of the integrand at X equals 0 and 1.0. 


Substitution of equation (Bl) into equation (36a) results in 
(after integration) 



Equation (B2) gives X n <_ as a function of the parameters I x and £. 

* * e 

The equation has been solved for a wide range of I x and \ , and the 
results are plotted in figure 16. 2 


* . 

For infinite, blade-element solidity, Ii_ and £ become 1.0 
so that (from fig. 16) X^ 5 becomes 0.5 (the value assumed, in effect, 

for the first approximation) . As the solidity decreases from infinity, 
I x and | became different from 1.0 and X Q 5 becomes different from 

0.5. In the second approximation of the blade-element-design method, 
the value £ can be estimated but the value of I* cannot be deter- 

a 

mined conveniently. However, figure 16 shows that for all values of 
| (but especially for values of \ between l /2 and 2 ) the value of 
Xq 5 given by Ii_ equal to 1.0 (that is, assuming a linear variation 

in I with X) il sufficiently accurate for the approximate methods of 
this report, provided ijf lies between 0.8 and 1.2 (that is, provided 

the integrand deviates oMly moderately from a linear variation, as 
assumed in the second approximation, for high-solidity blade elements) . 
And in any event, the values of Xq 5 s0 determined are more accurate 

than the value of 0.5 assumed by the first approximation. 


For the two numerical examples of the design method in this report, 

0.98 <lj£ <1.12 and 0.8 <| <2.1. These ranges in I* and | resulted 

a F 

in values of Xq 5 between 0.47 and 0.59. Such a range of Xq 5 is 

typical of high- solidity blade-element profiles with minimum radii of 
curvature at least approximately equal to the channel width between 
profiles, (in reference 13, for example, the values of Xq 5 vary 



NACA TN 2408 


51 


between 0.50 and 0.58.) Therefore, for high-solidity blade- element 
profiles -with relatively large radii of curvature that result in 

integrands Q cos that deviate only moderately frcm linear 

variations in the X- direction, the value of Xq 5 depends more on the 
magnitude of the integrand at X equals 0 and 1.0 (measured by |, 
fig. 16), than on the magnitude of the deviation (measured by I x , 

\ I* 

fig. 16). Thus, the value of Xq .5 determined by assuming a linear 

variation in I (that is, -£- Q cos 3*) with X is sufficiently 
x Po / 

accurate for the purposes of this report. 
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APPENDIX C 


VELOCITY DISTRIBUTION IN CIRCUMFERENTIAL DIRECTION ACROSS 


CHANNEL BETWEEN BLADE ELEMENTS 


The variation in the relative velocity Q in the circumferential 
direction across the c hann el Between high-solidity Blade elements is 
determined from considerations of the absolute irrotational motion of 
the fluid. Consider a fluid particle on a surface of revolution 
(fig. 5) . In the absence of entropy gradients the absolute motion of 
the fluid particle is irrotational so that 

0 = 35 + Qe) R ] “ is [%] 

which, from equations (ll) , (12), and (30), becomes 

2R*Urn sin a Qq sin a * dQg a Bq M 

°" — sj — + — ^ +K s 


or, from equations (6a) and (6c), and with the assumption that for 
high-solidity blade elements the variation in p’ across the channel 
is Kmfll ~i anfl that this variation can be neglected in the differential 
equation so that p ' Z p^ is a function of M only 5 


0 « 


2R*Ui 


'T 




sin a Q sin pA c sin a * 

+ £2 + R*Q cos p: 


dp 


0.5 


Bt 


0.5 dM 


T,* ' O' dQ 

R sxn p 0>5 ^ 


CT COS 


Po. 


Bq 

BX 


Along the strea m l in e , 

aq _ Bq + 

aM ii dM sm ax 40 ^ 
and from equations (2a), (2b), (ll) , (12) , and (30), 

Bq dQ _ g tan Pq. 5 Bq 
<5M dM “ 5 ' R* BX 


(Cl) 


(C2) 
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The derivative dQ/dM varies in the X-direction across the 
channel between "blade elements. The values dC^/dM and dQ-^/dM 
at X equals 0 and 1.0, respectively, are known from the prescribed 
variations in Q, q and with M. For high- solidity blade elements', 
it is assumed that the variation in dQ,/dM with X can be approximated 
by a straight line so that 


where 


dQ 

dM 


E + F X 


E = 


F = 


dQ 0 

dM 

<% 

dM 


dQo 

dM 


(C3) 

(C3a) 

(C3b) 


The differential equation (Cl) when combined with equations (02) and 
(C3) becomes 


|| - AQ - B(E+FX) 


C = 0 


(04) 


A = 


5' 


„ cos p' _ (sin pi c B ~~ + R* cos p’ _ — ■§— 
a K 0.5 \ ^0.5 R^, M 0.5 dM 

B = ~ R* cos P r sin p* 

O u.5 0.5 




(05) 

(06) 
(07) 


Equation (C4) can be integrated in the X-direction to give 



where the constant of integration has been dete rmin ed so that 
Q equals Qq at X equals 0. The distribution of Q circumfer- 
entially across the channel between blade elements is given by 
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equation (C8) . The coefficients B, C, ~E, and F are given by equa- 
tions (C6 ) , (C7 ) , (C3a) , and (C3b) , respectively. The coefficient A 
is selected so that Q is equal to Q-^ when X = 1.0. Thus from 
equation (C8), 


h ~ 

Qq exp 

(A)] + C 

i • 

- exp (A) 

A J 

i 

(-1) 

rl - exp i 

L A 

[H 



(C9) 


which is solved hy trial and error to obtain the proper values of k, 
that is., the value of A that satisfies the known value of B. 


For small values of A, equations (C8) and (C9) become 

f 3 

Q = exp (AX) + BEX + C + (BEA + BF + CA) ( |y + ~ + + 


and 


B 


" Qq exp (A)] 


C " 04 l 21 + 31 + 41 + 


•) 


E + (ea + f) + jr + jfc + • • •) 


(C8a) 


(C9a) 


which forms for equations (C8) and (C9) eliminate "small differences 
of large numbers" that otherwise appear in the numerical calculations 
when A is small ( |A| <0.2) . 
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APEEUDIX D 


consiMT Ap in direction of 0 


It is proposed to show that Ap is approximately constant in the 
direction of 0 (that is, in the direction of Z) if (tan p - tan pQg) 
is small in the direction of 0 (a condition that is approached in high- 
solidity blades with gradual variations in blade-element-profile thick- 
ness) and if Ap is small. Consider the magnitude of Ap at some 
point (X,M) if, for the same value of M, the initial flow direction 
Pq g at Xq g is changed an amount APq with Xq g and 6 remain- 
ing unchanged throughout the flow field. 

From equations (ll), (12), and (30), 


0 - 


'0 



(Dla) 


and in particular 


S 0.5 " 0 O ~ *0.5 


V 


From equations (Dla) and (Dlb), 


0 - Vs - V (x 6 - *0.5 B > 


(.mb) 


which, after differentiating with respect to M and multiplying by R, 
becomes 

tan ft - tan = 2* fel - d(X °-5 &) ) (B2a) 

K h 0.5 a \ <3M dM / v 

Also, for the initial values of p' with the same values of X and 6, 
tan - tan ^ 

Because X and 6 remain the same, equation (D2b) is subtracted from 
equation (D2a) to give 

tan p - tan p* = tan p Q .g - tan p^ (D3) 
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But from equation (39) , 


= P " Ap 


so that equation (D3) Becomes 

tan p - tan (p - Ap) = tan P 0>5 - tan (P Q>5 - Ap Q>5 ) 

which for small values of APq ^ Becomes 

tan p - tanAp „ „ tan P 0.5 " ^0.5 

tan P “ x + t an p tan Ap ~ tan Pq. 5 “ 1 + Ap 0<5 tan p Q>5 


(D4) 


But 


tan P Q = tan p - (tan p - tan p n ,_) 


0.5' 


and, neglecting powers and products of (tan p - tan P Q ) and AP Q 
which are assumed small, equation (D4) Becomes 

tan p - tan Ap _ tan p - Ap Q5 

1 + tan p tan Ap ~ 1 + Ap tan p 

0.5 

from which 

tan Ap ~ Ap 0 5 
or. Because Ap^ ^ is small, 

Ap « Ap Q>5 (D5) 
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TABLE I - FIRST NUMERICAL EXAMPLE 


[Plane, two-dimensional cascade: |3y, 0°j Py, 60° j Qy, 0.750} 

Qy, 0.290} Cp, 0.278} 0 , 1.215]] 


M 

^av 

AQ 

Qo 

% 

AP 

(Q ^av 

P 

5 

P 0.5 

6 

X 0.5 

First approx- 
imation 

Second approxima- 
tion 

0 

0.348 

0.200 

0.248 

0.448 


0 

0 

0.850 

2.9 

0.890 

0.559 

.1 

.396 

.279 

.256 

.535 

-.141 

.036 

5.3 

.765 

7.2 

.816 

.580 

.2 

.448 

.338 

.279 

.617 

-.189 

.088 

11.3 

.701 

12.3 

.751 

.586 

.3 

.503 

.378 

.315 

.692 

-.230 

.153 

17.7 

.661 

17.7 

.702 

.581 

.4 

.559 

.398 

.360 

.757 

-.259 

.230 

24.3 

.642 

23.5 

.673 

.568 

.5 

.611 

.398 

.412 

.810 

-.275 

.314 

30.9 

.644 

29.6 

.666 

.549 

.6 

.659 

.378 

.470 

.848 

-.273 

.400 

37.3 

.667 

35.9 

.680 

.527 

.7 

.699 

.338 

.530 

.868 

-.252 

.482 

43.6 

.711 

42.1 

.716 

.506 

.8 

.730 

.279 

.590 

.869 

-.212 

.555 

49.5 

.779 

47,9 

.769 

.490 

.9 

.748 

.199 

.648 

.847 

-.154 

.613 

55.0 

.874 

52.8 

.834 

.480 

1.0 

.751 

.100 

.701 

.801 

-.078 

.650 

59.9 

.996 

56.4 

.901 

.474 
























TABLE II - SECOND NUMERICAL EXAMPLE 

jjtLxed-flow impeller far centrifugal compressor: X^, 0.2j X-jy 1.939; cp, 1.05; cr, 2.56l] 


L* 

M 

R* 

H* 

a 

Qav 

Aft 

Qo 

Qi 

AP 

Wav 

1 

P 

& 

^0.5 

B 

*0.5 

First ap- 
prrnH-mfit-1 nn 

Second approxima- 
tion 


0 

0 

0.5 

8.0 

0 

0.500 

0.100 

0.450 

0.550 

-0.060 

-0.350 

-44.4 

0.853 

-44.5 

0.860 

0.497 

.1 

.078 

.503 

7.955 

3.5 

.502 

.111 

.447 

.558 

-.067 

-.289 

-35.2 

.739 

-34.9 

.741 

.506 

.2 

.158 

.513 

7.817 

7.1 

.509 

.143 

.438 

.581 

-.089 

-.230 

-26.9 

.661 

-26.5 

.664 

.517 

.3 

.239 

.530 

7.580 

10.8 

.521 

.195 

.423 

.618 

-.125 

-.159 

-17.8 

.596 

-17.4 

.603 

.529 

.4 

.324 

.555 

7.231 

14.6 

.538 

.263 

.406 

.670 

-.178 

-.065 

-6.9 

.542 

-6.5 

.555 

.542 

.5 

.413 

.589 

6.753 

18.6 

.561 

.342 

.391 

.732 

-.247 

.063 

6.4 

.507 

6.8 

.529 

.553 

.6 

.509 

.634 

6.118 

22.9 

.592 

.420 

.382 

.803 

-.334 

.226 

22.4 

.505 

22,9 

.537 

.555 

.70 

.612 

.694 

5.290 

27.6 

.634 

.482 

.393 

.875 

-.433 

.414 

40.8 

.573 

41.6 

.617 

.541 

.75 

.668 

.729 

4.787 

30.1 

.660 

.498 

.411 

.909 

-.484 

.501 

49.4 

.647 

50,4 

.698 

.528 

.800 

.726 

.770 

4.215 

32.7 

.690 

.498 

.442 

.939 

-.531 

.568 

58.4 

.788 

59.6 

.843 

.510 

.825 

.757 

.793 

3.900 

34.1 

.708 

.490 

.463 

.953 

-.551 

.622 

61.5 

.863 

62.7 

.910 

.525 

.850 

.788 

.817 

3.563 

35.5 

.727 

.476 

.489 

.965 

-.566 

.647 

62.8 

.906 

63.8 

.940 

.547 

.875 

.821 

.843 

3.204 

36.9 

.749 

.455 

.521 

.976 

-.575 

.658 

61.5 

.884 

61.9 

.921 

.581 

.900 

.854 

.870 

2.819 

38.4 

.772 

.425 

.560 

.985 

-.574 

.653 

57.7 

.815 

56.5 

.825 

.570 

.925 

.889 

.899 

2.409 

40.0 

.799 

.386 

.606 

.991 

-.560 

.631 

52.2 

.751 

49.8 

.745 

.553 

.950 

.925 

.931 

1.971 

41.6 

.828 

.336 

.660 

.996 

-.526 

.588 

45.2 

.719 

41.0 

.704 

.520 

.975 

.962 

.964 

1.502 

43.3 

.862 

.275 

.725 

.999 

- .466 

.523 

37.4 

.747 

32.6 

.744 

.463 

.9875 

.981 

.982 

1.255 

44.1 

.880 

.239 

.760 

1.000 

-.423 

.482 

33.2 

.800 

29.8 

.829 

.430 

1.0 

1.0 

1.0 

1.0 

45.0 

.900 

.200 

.800 

1.0 

-.370 

.439 

29.2 

.906 

25.8 

.958 

.421 
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Figure 3. - Surface of revolution vith coordinates and velocity c components. 
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Center line (between 
adjacent streamlines 
in meridional plane) 
used to generate 
surface of revolution 


Rd0 
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Blade -element spacing 


at tip, %(A8) = — 


t -boundary 

, (M = 1.0) 


Me / 9 


Blade-element / (/ 
camber line^/ / 



Blade - element spacing 
- / R(AS) = 


f(l-6)R* / Channel vldth 


RCOi-Oo) = R* I 


n-boundary 
(M = 0) 


Figure 6. - View of passage between two blade-element profiles on 0M-plane 
showing parameters used in first approximation. 
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t-boundary 



N 


Figure 11. - Comparison of blade-element profiles obtained by first and second approximations for 

first numerical example. 
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Figure 12. - Circular arc used to generate surface of revolution for second 
numerical example. Variation in H* along circular arc also shown. 
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(a) Surface of revolution. 

figure 13. - Blade-element profile for second numerical example. 
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1.0 


.8 — 


• 6r~ 


L* 


• 4h 


• 2 — 


0L 

-.4 


x Stream function, i|r 

t -boundary 0 .2 .4 .6 .8 1.0 



(b) Transformed L*8*-plane shoving streamlines 
used in stream-filament check. 


Figure 13. - Concluded. Blade-element profile for second 

numerical example. 
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i l i i i i i i i i i 

0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 


L* 

Figure 14* - Comparison of prescribed velocities for second numerical example with velocities obtained 
by stream- filament check' of profile obtained by second approximation* 
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Figure 16. - Effect of ij (measure of deviation of parabolic variation in X 


f> * * 

from linear variation) and | (ratio of 1^ to I q ) upon ratio Xq.5» 


NACA - Langley Field, V*. 



